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The stellar populations of massive galaxies 
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Abstract. I present a brief review of the stellar population properties of massive galaxies, 
focusing on early-type galaxies in particular, with emphasis on recent results from the ATLAS 3D 
Survey. I discuss the occurrence of young stellar ages, cold gas, and ongoing star formation in 
early-type galaxies, the presence of which gives important clues to the evolutionary path of 
these galaxies. Consideration of empirical star formation histories gives a meaningful picture 
of galaxy stellar population properties, and allows accurate comparison of mass estimates from 
populations and dynamics. This has recently provided strong evidence of a non-universal IMF, 
as supported by other recent evidences. Spatially-resolved studies of stellar populations are 
also crucial to connect distinct components within galaxies to spatial structures seen in other 
wavelengths or parameters. Stellar populations in the faint outer envelopes of early- type galaxies 
are a formidable frontier for observers, but promise to put constraints on the ratio of accreted 
stellar mass versus that formed 'in situ' - a key feature of recent galaxy formation models. 
Galaxy environment appears to play a key role in controlling the stellar population properties 
of low mass galaxies. Simulations remind us, however, that current day galaxies are the product 
of a complex assembly and environment history, which gives rise to the trends we see. This has 
strong implications for our interpretation of environmental trends. 
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1. Introduction 

Much of what has been written about the stellar population content of nearby massive 
galaxies has focussed on that of early-type galaxies specifically, by obvious virtue of their 
observational simplicity, having relatively high surface brightnesses, simple morphologies, 
and very little dust or emission to complicate the interpretation of their spectra. I will 
not deviate from this well-worn path in this review, but point out that 'massive galaxies' 
need not be synonymous with early-type galaxies alone, which only dominate the stellar 
mass function above 3.el0 Mq (e.g. Baldry et al. 2006). That said, a large portion of the 
stars in the local Universe do exist in early-type galaxies, mostly in the form of old stellar 
populations that carry the fossil evidence of earlier evolutionary epochs. In this review, 
I focus on the properties of early-type galaxies, with a bias towards recently published 
results from the ATLAS 3D Survejjfj (Cappellari et al., 2011). 

2. Colour bimodality 

The colour-magnitude diagram of galaxies (in particular that from SDSS, e.g. Kauff- 
man et al., 2003, Baldry et al., 2004) is an excellent tool to summarize our current knowl- 
edge and ideas of galaxy evolution, being a snapshot of the cumulative star formation 
history (color) and mass assembly (magnitude) for the population. In an evolutionary 
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context, it is inevitable to ask how galaxies may move around in this parameter space as 
a function of time, and what processes give rise to the main features of the distribution. 

Within the hierarchical framework of ACDM, we assume that galaxies largely follow 
their host dark matter haloes, and increase in mass over time. Likewise, the stellar popu- 
lations will evolve over time, becoming redder. New stars may form if sufficient fuel exists 
and conditions allow it, producing or preserving bluer colors. Various tracks through the 
color- magnitude distribution have been proposed, most of which consider that blue galax- 
ies are the building blocks for most red-sequence galaxies, with evolution along the red 
sequence also necessary to produce the most massive galaxies, which have insufficient 
blue counterparts of similar mass. The gradual reddening process of gas exhaustion and 
stellar evolution, however, do not obviously give rise to the observed bimodal color distri- 
bution. Creating an accelerated transition from blue to red colors can be achieved with 
the addition of a mechanism which 'quenches' star formation (Faber et al., 2007). 

3. Feedback 

One suggested mechanism for quenching star formation is that of energy input into 
the star-formation gas reservoir. The form that this energy takes is not currently clear, 
nor is the precise source of energy input (AGN or stellar winds). Circumstantial evidence 
exists, however, to support the cessation of star formation due to the onset of black hole 
accretion. Specifically, the work of Schawinski et al. (2007) found that the most active 
AGN galaxies lie in the intermediate region between red and blue galaxies but have only 
small amounts of recently-formed stars, strongly suggesting that the AGN activity is 
associated with the quenching of star formation. 

An example of a direct connection between nuclear accretion and the removal of star 
formation fuel comes from the recent discovery of a massive molecular gas outflow from 
the early-type galaxy NGC 1266 (Alatalo et al., 2011, Davis et al., 2012). This otherwise 
typical lenticular galaxy shows molecular gas flowing away from the central regions faster 
than the escape velocity, losing gas faster through the outflow than via the apparent star 
formation rate. X-ray and radio signatures suggest the presence of an AGN, which would 
be an obvious power source. This object appears somewhat unique within the volume- 
limited parent sample of 260 early-type galaxies it was taken from, observed as part of 
the ATLAS 3D Survey. However, it is a bright and easily detected case, and less obvious 
examples may be lurking below the sensitivity limits of current surveys. 

4. Cold gas in early-type galaxies 

The largest and most sensitive surveys targeting early-types in CO (single dish, J=l-0, 
1-2) and HI (interferometric, 21 cm) were recently completed as part of the ATLAS 3D 
Survey (Young et al., 2011, Serra et al., 2012), building on previous investigations (Mor- 
ganti et al., 2006, Young et al., 2008). By virtue of being a complete sample, ATLAS 3D 
has provided secure detection fractions for HI and CO down to sensitive mass limits 
(around 10 7 M ). Both tracers show numerous cases where the gas mass and gas frac- 
tion are comparable to that of spiral galaxies. More than 20% of early-type galaxies show 
detectable cold-gas content in either tracer, with a strong correspondence of HI and CO 
in the same objects, and little dependence on galaxy mass (see Figure 1). 

CO interferometry was obtained for the 30 brightest detections in the ATLAS 3D Survey 
(Alatalo et al., 2012), allowing the spatial distribution of the CO gas to be investigated. 
The extent and star-gas kinematic alignment derived from these data have revealed inter- 
esting constraints on the origin of the molecular gas in these objects, implying important 
differences in the molecular gas properties in different environments (Davis et al., 2011, 
see also Davis et al. this volume). Specifically, objects in the Virgo cluster (the densest 
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Figure 1. Left: Molecular gas masses for the ATLAS 3D sample as a function of dynamical mass, 
from Young et al. (2011). Circles and crosses are detections and non-detections respectively. 
Righ t: Neutral gas mas s for ATLAS 3D galaxies as a function of absolute K-band magnitude, 
from Serra et al. (2012) Filled and open symbols are detections and non-detections respectively. 



environment of the survey) almost exclusively show aligned, co-rotating CO gas, unlike 
the field which shows a similar occurrence of aligned and misaligned gas kinematics. 
This supports a picture where field galaxies continue to accrete gas via minor mergers 
with in falling systems, giving a range of resulting orientations, whereas in the cluster 
environment, this channel is apparently suppressed, and an internal origin seems likely. 

Detection rates of CO in Virgo are also as high as in the field (Young et al., 2011). This 
contrasts with the HI detection rate, which is almost nil within the cluster, suggesting 
that the denser CO gas is more robust to the tidal interactions and ram pressure of 
the cluster environment than the HI. HI morphology also shows a great variety, with 
disks, resonant rings, and tidally-disturbed structures, all of which can exist on 10s of 
kpc scales (Serra et al., 2012). The CO, by contrast, is typically found in regular rotating 
disk structures within the central few kpc (Alatalo et al., 2012, Davis et al., 2012). 

5. Recent and ongoing star formation in red-sequence galaxies 

The presence of neutral gas, and in particular molecular gas, at levels of up to 10% by 
mass is a strong indication that these objects are currently experiencing star formation. 
This is consistent with findings from the SAUR0N Project, which looked for excess flux in 
the 8/im Spitzer/IRAC band, indicative of PAH emission from star formation (Shapiro et 
al., 2010). The rates of star formation found are typically low (around 0.1 M /yr) with 
some objects getting closer to spiral-like values (1 M Q /yr). The star formation rate does 
not depend strongly on galaxy mass, although the morphology of the star-forming region 
does. More massive galaxies generally show star formation within disks that have small 
scale lengths compared to the half-light radius of the stars, whereas going to lower mass 
galaxies, the star formation region can be as extended as the stellar light. This creates 
very different impacts on the color and, more sensitively, the mean luminosity-weighted 
age of the galaxy. 

Evidence of young stars in early-type galaxies comes from the first detailed studies of 
stellar populations using absorption lines and population synthesis models (e.g. Gonzalez, 
1993; Trager et al., 2000). Young mean ages inferred from integrated optical light is 
straightforwardly interpreted as a 'frosting' of young stars superimposed on a body of 
mostly old stars. Later studies refined the picture to show, for large samples, that young 
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galaxy ages were more common in objects with lower velocity dispersion (Thomas et al., 
2005), and indeed the two parameters are tightly correlated (e.g. Graves et al., 2009). 
The presence of young stars is most easily seen at bluer wavelengths, making the ultra 



violet (UV) a particularly sensitive tracer of stars formed in the past Gyr. Kaviraj et 



al. (2007) showed the dramatically increased dispersion of the red sequence when using 
UV-optical colors, driven primarily by the presence of young (<J 1 Gyr) stars in more 
than 30% of the present day early-type galaxy population. 

6. Star formation histories of early-type galaxies 

As high-quality stellar libraries and associated stellar population models have become 
available (e.g. Sanchez-Blazquez et al., 2006; Vazdekis et al. 2012), it is now possible to fit 
sizable optical-NIR wavelength ranges at moderate spectral resolution with spectra that 
each represent a population of fixed age and metallicity, with additional free parameters 
like the assumed initial mass function and elemental abundances (Conroy & van Dokkum, 
2012). This allows the full information content of the spectrum to be used, and avoids 
changing the observed data to accommodate the specific model framework by instead 
adapting the models to the data. 

A number of approaches have been employed to use such models and spectral fitting to 
derive the star formation histories of galaxies, typically either selecting a fixed functional 
form (such as an exponential), and solving for the function's free parameters (e.g. Pforr 
et al. 2012); or selecting the best fitting scenario from a large library of different trial 
histories, which can take arbitrary forms (Aquaviva et al., 2011). Another approach is 
to take a regularized linear combination of model components and find the 'smoothest' 
distribution of free parameters that can still adequately fit the data (Ocvirk et al., 2006). 

7. Galaxy mass-to-light ratios 

As well as being instructive from a galaxy evolution perspective, knowing the distri- 
bution of population parameters within a galaxy can have a large impact in how we 
interpret the stellar light in terms of mass, via the stellar mass to light ratio. For exam- 
ple, a galaxy may have absorption line strengths that are fully consistent with a single 
stellar population (SSP) age of 3 Gyr. The same combination of indices may be equally 
well fitted by combination of two different SSPs, say 95% by mass of 13 Gyr and 5% of 
1 Gyr. These two scenarios would differ in inferred mass by almost a factor of 3. 



This effect was shown in Cappellari et al. (2006) where the population and dynamical 



M/L properties of 24 early-type galaxies, all with high-quality integral- field spectroscopy 
from the SAUR0N Survey, were compared. Younger galaxies deviated systematically from 
the identity line towards lower population-based M/L. By modeling these younger galax- 
ies using two population components, one older (12 Gyr, 90% by mass), and one younger, 
the dynamical- and population-based M/L values came into agreement. 

An analogous study was recently made by Cappellari et al. (2012a, 2012b) using the 
ATLAS 3D Survey of 260 early-type galaxies, again modeling the stellar dynamics and 
populations derived from integral-field spectroscopy covering to around one effective ra- 
dius for each galaxy. Star-formation histories were derived for each galaxy using the 
regularized spectral fitting method of Cappellari fc Emsellem (2004) | fitting for the age 



and metallicity distribution in each object. From these distributions of mass contribu- 
tions in each age and metallicity bin, it is possible to compute the mass in stars, as well 
as the light output, giving effectively the mass-weighted M/L. 

For the dynamical M/L estimates, it is first necessary to describe the stellar light 
distribution. It is important to do this as accurately as possible, as it will determine 
the assumed stellar mass distribution, which will in turn affect the spatially resolved 
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Figure 2. Ratio of stellar M/L derived from stellar dynamical modeling and stellar population 
models, plotted as a function of dynamical M/L. Symbol colours correspond to velocity dis- 
persion. Horizontal lines indicate where the points would lie for different canonical IMFs. The 
six different panels illustrate the best fitting results under different dark matter assumptions, 
as indicated in each panel. In all cases, the ratio of the two M/L estimates is not constant, 
and shows the same systematic deviation towards higher mass normalisations in higher mass 
galaxies. From Cappellari et al. (2012a). 



kinematics predicted for a given M/L. Early- type galaxies show a wealth of detailed 
structure in their spatially resolved surface brightness distributions, even though they 
are often assumed to follow computationally convenient, spherically-symmetric profiles. 
Instead, the 260 SDSS r-band galaxy images were fitted in detail using the Multi-Gaussian 
Expansion (MGE) approach (Monet et al., 1992; Emsellem et al., 1994). Even imposing 
the constraint of axisymmetry, the light distribution is fitted to high precision, matching 
substructures such as disks, boxiness, and variations in radial profile shapes that are 
missed in simpler models. 

The dynamical M/L estimates were computed using the Jeans Anisotropic MGE 
(JAM) approach of Cappellari (2008) | which exte nds the semi-isotropic (/ = f(E,L z ), 
two integral) formalism of |Emsellem et al. (1994) to the case where the velocity ellipsoid 
is aligned with the cylindrical coordinates (R, z), and the anisotropy is parameterised as 
/3 Z = 1 — v\lv\. The free parameters are the M/L, j3 z , and the inclination of the system, 
i. Under these assumptions, the Jeans equations are solved to produce an estimate of the 
projected second velocity moment V rms = s/v 2 + <r 2 as a function of position, which is 
compared to the observed kinematics, and the best-fit case found. The resulting fits to 
the high-quality kinematic data tightly constrain the M/L. 
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8. Measuring variations in the IMF 

In addition to the effects of young stellar populations on the inferred M/L, jCappel- 



lari et al. (2006) also showed that, adopting a universal form of the stellar initial mass 
function (IMF), it is not possible to have the M/L from populations and dynamics agree 
for all galaxies. Specifically, there was a weak trend that higher-mass galaxies showed 
increasing dynamical M/L with no corresponding increase in population M/L - i.e. the 
stellar kinematics showed a mass component in these galaxies that is unaccounted for 
in the light. Two possibilities are: (1) that the dark matter fraction is different in the 
central parts of these galaxies; or (2) the initial mass function of the stellar populations 
is different. However, the effects of these two possibilities are degenerate for the data and 
modeling employed. Interestingly, this same finding was later confirmed using different 
data and an independent approach, via gravitational lensing (Treu et al., 2010). 

With ATLAS 3D , the sample size increases ten- fold over the SAURON study, with compa- 
rable 2D kinematic and imaging data. The JAM modeling approach is computationally 
much more efficient compared with more general methods , and is readily applied to large 
samples. Integral field data are also essential here, as the inclination and anisotropy vary 
the shape of the V rms distribution away from the principle axes, which would not be seen 
with major- or minor-axis long slit data. 

The fit to the stellar light and 2D V rms is sufficiently constraining that, for the first 
time, it is possible to break the degeneracy between the stellar mass distribution and the 
dark matter component, thus allowing trends in these properties to be tracked across the 
galaxy sample. This was done by assuming several different models for the distribution of 
the dark matter, each motivated by current theoretical models. Figure [2] shows the results. 
If the estimates from both methods were consistent, the points would all lie horizontally. 
Since the dark matter variations are accounted for, this variation must come from the 
initial mass function of the stellar population (which is held fixed in this diagram). 

This finding supports recent claims (van Dokkum & Conroy, 2010, Conroy & van 
Dokkum 2012) that show spectral evidence for large numbers of faint, low-mass dwarf 
stars in the centers of early-type galaxies. This claim has been generally supported by 
numerous other similar studies (e.g. Spiniello et al., 2012; Ferreras et al., 2012; Smith et 
al., 2012b). We note that these spectral studies are required to distinguish the slope of 
the IMF, as gravitational studies alone cannot distinguish between dark mass from faint 
low-mass stars, or remnants from massive star evolution. 

9. Spatially-resolved stellar populations 

The variations of stellar populations within massive galaxies indicate the diversity 
of processes in their evolution, and are consistent with a number of possible formation 
scenarios. Once again, the power of integral- field spectroscopy gives a unique view of these 
properties, allowing spatial structures in the populations to be unambiguously associated 
to other features, such as kinematic components, dust lanes, gas disks, etc. In this way, 



Kuntschner et al. (2010) showed that the presence of young stars is often accompanied 
by the presence of a rapidly-rotating disk component in the kinematics. In some cases, 
these components also show disks of molecular gas that are co-spatial with the young 
stellar disk (Young et al., 2008), showing that gas has settled onto the equatorial plane 
of the system and is forming new stars. 

Metallicity gradients are essentially ubiquitous in early-type galaxies of all masses, and 
occur in the sense that metallicity decreases in the outskirts, falling approximately as 
a power law with radius. A more physical interpretation of these gradients comes from 
considering that the metallicity is tracing a change in the gravitational potential, which 



can be characterized by the escape velocity. A recent study by Scott et al. (2009) has 
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Figure 3. Left: Radial line-strength gradients for NGC 3379, fr om |Weij mans et al. (2009) 

389, from[Coccato et al. (2010) ~ 



Right: Radial profile of SSP parameters for NGC 4889, 



shown that early-type galaxies collectively follow trends in there local stellar population 
- escape velocity relations. This is suggestive that the formation of metals does not 
care about the specific galaxy as much as what the local potential is. In a situation 
where stellar winds control the star-formation process, deeper potentials could retain 
more enriched gas, leading to a strong connection between the metallicity and the local 
potential. Dry galaxy merging would disrupt this connection, and so the tightness of the 
local-global trends limits the number of major dry mergers to approximately 1.5 over a 
Hubble time (Scott et al. 2009). 

A somewhat different scenario for building metallicity gradients in was proposed by 
Oser et al. (2010)] based on simulations. Here the central regions of the most massive 
galaxies form very early on, through rapid gas inflow and subsequent collapse, quickly 
forming a moderately massive, compact object. These then proceeded to accrete lower 
mass galaxies. The merging process is such that most of the stars of an accreted satellite 
infall only as far as where the satellite's binding energy equates that of the parent poten- 
tial. Thus stars from lower-mass galaxies (which in turn have lower metallicity due to the 
effect of stellar winds) remain at larger radii, naturally creating a metallicity gradient. 

Studies of metallicity profiles at very large radii are difficult, as the galaxy surface 
brightness becomes very low - orders of magnitude below the sky level. Nevertheless, 
innovative observations, such as deep long-slit integrations (Coccato et al., 2010), use of 
integral-field spectroscopy as a large single aperture (Weijmans et al., 2009, Greene et 
al., 2012), or using multi-object spectroscopy to collect light over a 2D region (Foster et 
al., 2009) or from proxy tracers like globular clusters (Arnold et al., 2011), can measure 
line-strengths out to 3-4 effective radii (see Figure 3) . In general, changes from the central 
gradient are subtle, and current small samples preclude general conclusions. 

Colour profiles offer a promising technique to study these extremely faint outer regions. 
Deep imaging of massive galaxies often reveals direct signs of interactions (Due et al., 
2011). Stacking carefully selected shallower data can also probe depths out to 8 effective 
radii (La Barbera et al., 2012) - much further out than is possible with current spec- 
troscopic instruments. In addition to the difficulties of accurate sky subtraction and flat 
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Figure 4. L eft: Stellar age versus mass and projected radial distance from the centre of the Coma 
cluster, from |Smith et al. (2012a)l Right: Fraction of red galaxies versus mass and environmental 
density, from ]Fenget^L(27n0)p rhe right-hand figure has been inverted to show density varying 
in the same direction in both figures. The matching thick hand-drawn line on both figures 
indicates the similar behavior of age and red fraction with mass and environmental density. 



fielding at these very low surface brightnesses, the effects of chromatic PSF differences 
is very important, as the faint halo of the PSF spreads low-level light across large scales, 
and varies with each filter. Secondly, the stellar population models suffer degeneracy be- 
tween age and metallicity in broad-band photometry, and generally lack sensitivity to 
subtle metallicity changes. However, the extensive baseline in radius, and rich spatial 
information make deep imaging an important complement to spectroscopy. 



10. Environment 

Much has been written on the role of environment in influencing galaxy evolution, and, 
in particular, the star-formation process within galaxies. Rich clusters are the largest over- 
densities of baryonic matter in the Universe, and in a hierarchical framework, assembled 
mass more quickly than their surroundings. But does the act of being in a cluster impose 
intrinsic stellar population differences on member galaxies compared to the field? 

The influential study of Thomas et al. (2005)| found a difference in the mean ages of 



massive galaxies in high-density and low-density regions independent of galaxy velocity 
dispersion. In a more recent study, using similar methods for a much larger sample 



of galaxies but with lower quality data, Thomas et al. (2010) revised this picture to one 



where the more massive population has stellar ages that are independent of environment. 
Instead, the difference is attributed to the significance of the low mass 'rejuvenated' 
population, which are more common in low-density environments. 



A recent study by Smith et al. (2012a) addressed the question of environment using 
SDSS spectroscopy of the Coma cluster. This study found that, at high masses, the 
galaxy SSP-equivalent ages are only weakly sensitive to their position in the cluster, in 
the sense that all the most massive objects are old regardless of their position in the 
cluster. By contrast, the age of low-mass galaxies vary strongly with radius from the 
cluster centre. The combined effect is a reduction in the mean stellar age of galaxies with 
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increasing distance from the cluster centre, which is mostly driven by the dependence of 
the lower-mass objects to the cluster-centric radius. 

We show in Figure 4 that this picture is similar in essence to that of the apparent 
'separability' of mass and environment with respect to the quenching of star formation, 



proposed by Peng et al. (2010) Whereas the specific mechanism for 'mass quenching' 
is rather unclear, the fact that environment has an effect preferentially on lower mass 
systems is expected, given the clear signs of interactions and stripping that characterize 
the experience of a gas-rich galaxy falling into a larger parent or cluster halo (e.g. Vollmer 

et al., 2010). 

As emphasized in Smith et al. (2012a) | the radial age dependence in Coma, driven by 
lower mass galaxies, can be understood when one considers that clusters are not well 
mixed and relaxed dynamical systems in the hierarchical paradigm. Indeed, for fixed 
halo mass, simulations show that the present-day position within the cluster is closely 
related (on average) to the infall epoch in the past (e.g. Smith et al., 2012a; De Lucia 
et al., 2012; McGee et al., 2012). Thus environmental quenching associated to a galaxy 
transitioning from a 'central' to a 'satellite' object naturally imposes a temporal and 
spatial signature on the stellar population properties of cluster galaxies. Mass assembly 



is closely connected to this process via accretion and merging, and De Lucia et al. (2012) 



therefore argue that the 'history bias' inherent in the observational snapshot we see at 
any epoch may render the question of mass versus environment ill-posed. 



11. Conclusions 

The main conclusions I wished to convey with this contribution are 

• The 'massive' galaxy population has a significant contribution from late-type spiral 
galaxies at all but the very highest masses. 

• Compelling circumstantial evidence connects AGN activity with the transition of 
galaxies from the star-forming blue cloud to the passive red sequence. However, detailed 
studies of objects 'caught in the act' of reducing their star-formation are needed to 
progress from correlation to causation. 

• Cold gas, young stars, and ongoing low-level star formation are frequently observed 
in early- type galaxies, and must be considered as part of their formation/evolution 

• The IMF varies systematically between low- and high- mass early- type galaxies. 

• Spatially-resolved observations are a unique way to accurately connect observational 
features across widely different datasets. Spatial distributions and kinematic orientations 
also offer unique constraints on the origins of galaxy components. 

• The outer envelopes of galaxies are new frontiers for testing galaxy formation models 

• Environment plays a key role in our picture of galaxy formation and evolution. 
Insights from galaxy formation models help us to see these trends as part of a complete 
formation path, which may not easily be inferred from instantaneous view we have. 
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